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Abstract 
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results are necessary for achieving low transverse momentum resummation at next-to-next-to- 
leading-logarithmic accuracy in this effective field theory approach with unintegrated distribution 
functions. 
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I. INTRODUCTION 



Factorization plays a crucial role in collider physics by allowing for a more predictive 
framework through the separation of perturbative and non-perturbative effects. For pro- 
cesses such as fully inclusive Drell-Yan production of lepton pairs, factorization expresses the 
theoretical prediction as the convolution of a perturbatively calculable partonic cross-section 
with universal non-perturbative parton distribution functions (PDFs). Large logarithms of 
the hard and soft scales are resummed by evolving the PDFs from Is.qcd up to the hard scale, 
of order the invariant mass of the final-state leptons, via the DGLAP renormalization-group 
equations. For more exclusive processes, restrictions on the final state can introduce new 
functions associated with intermediate momentum scales. In such cases, additional resum- 
mation is required and often new non-perturbative functions beyond the standard PDFs can 
arise. 

An important example of such an observable is the low transverse momentum {pt) dis- 
tribution of electroweak gauge bosons and the Higgs boson. It plays an important role in 
the precision measurement of the W^-boson mass, Higgs boson searches, tests of perturbative 
Quantum Chromodynamics (QCD), and in probing non-perturbative transverse momentum 
dynamics in the nucleon. For perturbative values of pt-, three distinct scales appear in this 
problem, M ^ pt ^ ^qcd, leading to a more intricate factorization formula. Large log- 
arithms of M/pt and Pt/ ^qcd must be resummed. In the non-perturbative region where 
Pt ~ ^QCDi new non-perturbative structures that probe transverse momentum dynamics in 
the nucleon appear. 

Low Pt distributions have been extensively studied in the traditional QCD literature 
|T2|, and have also been explored |]T3|, |T^ in the context of the Soft-Collinear Effective 
Theory (SCET) [|T3|-|T^. Recently, a new approach, based on SCET and on fully uninte- 
grated nucleon distribution functions, was developed in Refs. [^^-^, resulting in a new 
factorization and resummation theorem for low-p^ distributions. This approach allows one 
to predict the perturbative pt distribution entirely in terms of perturbatively calculable 
functions and the standard PDFs, avoiding the difficulties in matching the low and high 



Pt regions |Tl], ^ associated with treating the Landau pole in traditional approaches. The 

ch takes the schematic form 

H®Br,0Rn^S~\ (1) 



factorization theorem in this new approach takes the schematic form 



dpj. dY 

where H denotes a hard function, are Impact-parameter Beam Functions (iBFs), and 
S^^ is the Inverse Soft Function (iSF). All objects have well-defined operator definitions, 
as shown in Refs. [18, 19]. The hard function H is perturbatively calculable and encodes 
hard physics of the Drell-Yan production vertex. The iBFs are fully unintegrated nucleon 
distribution functions and the soft function S is given by the vacuum matrix element of soft 
Wilson lines. The iSF S^^ appears instead of S because of zero-bin subtractions p2|-p4 
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necessary to avoid double counting the soft region, as explained in Refs. [ITH], |T^. Similar 



soft-subtractions ||2^, appear in the formalism based on transverse-momentum dependent 



PDFs (TMDPDFs) p5|-p5|. For non-perturbative values of pr, the iBFs and the iSF are 
non-perturbative functions that encode the physics of non-perturbative pt emissions and 
transverse momentum dynamics in the initial state nucleons. For perturbative values oi px, 
the iBFs describe the evolution and shattering of the initial state nucleon into an initial 
state beam-jet of high energy pj- radiation. Analogous beam functions were first shown 



to arise in other contexts and correspond to a special case of the iBF. The iSF 

describes the emission of low energy p^ radiation from the initial states. In the standard 
approaches, rapidity divergences arise in perturbative computations of the TMDPDFs that 
are not regulated in dimensional regularization and are instead regulated with additional 
external regulators. In contrast, the iBF, which is more differential in momentum coordinates 
than the TMDPDF, can be computed in standard dimensional regularization with rapidity 
divergences regulated by the physical kinematics of the process. 

For perturbative values of pr, the iBFs can be perturbatively matched onto the standard 
PDFs, thus factorizing the non-perturbative dynamics of the initial state nucleon from the 
perturbative pt emissions. In this case, the factorization theorem for the cross-section, 
differential in the px and rapidity (Y) of the electroweak gauge boson, takes the form 

— — / dxi / dx2 



dprp dY Nc Jq Jq Jxj^ x[ J^^ x'2 

X i^|(xlX2Q^ /ig; [It) ^^"(a^i, X2, x\, x'^.Pt, Y, fiT)fr{x[,fiT)fs{x2, /ir), 

(2) 

where Q denotes the hadronic center of mass energy and the Transverse Momentum Function 
(TMF) function Q'^^'^ is given by 

g'^'''{Xi,X2,x[,X2,PT,Y,flT) = J J^^Jo[b±PT] J dt^dt^ In;qri^ ,tt , b±, fXr) 



X2 2^2^^ XiLJ 

(3) 

The functions Tn-qr and In-^qs are Wilson coefficients that arise from the perturbative match- 
ing of the iBFs onto the PDFs and are given by the finite part of the iBF computed in 
pure dimensional regularization. The leading-order (LO) and next-to-leading-order (NLO) 



expressions for the iBFs and the iSF were computed in Ref. |T9[ and were used to calculate 
the next-to-leading-log (NLL) perturbative pt spectrum for the Z-boson. A resummation 
at the next-to-next-leading-log (NNLL) level requires a computation of the iBFs and the 
iSF at next-next-to-leading-order (NNLO) in perturbation theory. 

In this paper, we take the first step towards achieving a NNLL resummation of the Drell- 
Yan pr-spectrum, using the effective field theory approach with unintegrated distribution 
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functions, by computing the soft function S that appears in Eqs. (||) and @ at NNLO. We 
perform several consistency checks on our calculation. Both the result for this exclusive soft 
function, and the techniques used in deriving it, should be of use in other investigations of 
resummation to high accuracy within effective field theory. Recently [P8|-^, two-loop results 
for a related soft function that appears in thrust distributions [|^-^ of e+e~ collisions were 
given, demonstrating the the arising need for studying higher-order corrections to multi-scale 
objects appearing in factorization theorems. Our paper is organized as follows. We formulate 
the problem and introduce our notation in Section We present both the techniques for 
and results of our calculation, including a comparison with known results, in Section |T|. 



Finally, we conclude in Section IV 



II. NOTATION AND NLO RESULTS 



The operator definition of the soft function in position space is given by 
S{b,f,) = ±TT{0\f[SlSn]{b) T[4^n](0)|0), 



(4) 



where Nc = 3 denotes the number of colors. Sn,n denote soft Wilson lines along the n^, n'^ 
directions respectively and are defined as 



Sn = Pexp 



ds n ■ Asix + sn) 



si = p 



exp 



ds n ■ Asix + sn) 



, (5) 



with analogous definitions for Sn and 5*1. The four-vectors n^,nf^ are light-like and satisfy 
n ■ n = 2. The symbols P, P denote path-ordering and anti-path-ordering respectively, and 
similarly T, T denote time-ordering and anti-time-ordering. 

One can define a hybrid soft function with light-cone momentum coordinates and position 
space impact-parameter coordinates as the Fourier transform of S{b, fi) with respect to the 
light-cone coordinates as 



S{q ,q^,b±,fi) 



db+db" 
167r2 



^iq-b+/2iq+b-/2 



S{b,fi), 



It is this hybrid soft function that appears in the factorization theorem in Eq. 
the full momentum space soft function is defined as 



(6) 

Similarly, 



5(g,/i) 



d^-^b 



(27r)^-2 

or equivalently, is related to the full position space soft function in the standard manner 

d'^b 



(7) 



S{q,n) 



(2vr)< 



(8) 



4 



We will present results for the exclusive NNLO position-space soft function S{b,fi) and 
the hybrid impact-parameter space soft function S{q~ ,q~^ ,b±) that appears directly in the 



factorization theorem of Ref. [0. The momentum-space soft function S{q) of Eq. is 
used in intermediate stages of the calculations, and is simple to derive using the presented 
formulae. We note that the position- space soft function of Eq. (|) by definition is equal to a 
gauge invariant soft function but evaluated in covariant gauges. In non-covariant or singular 
gauges such as the light-cone gauge, additional transverse gauge hnks are required in the 



definition of the soft function |S3[. These transverse gauge links are unity in covariant 



gauges and thus do not appear in the definition of Eq. (^). Similar arguments also appear 



in Ref. [^, |30| where the TMDPDFs were calculated in non-singular gauges. We restrict 



our analysis here to covariant gauges so that gauge invariance is fully respected. 
A. Renormalization 

Here we discuss the renormalization conventions for the hybrid impact-parameter space 
soft function S{q~ ,q~^ ,bj_, fi) that appears in the factorization theorem of Eq. (|]) and for 
the position-space soft function S{b,fi) of Eq. (^). We regulate infrared and ultraviolet 
divergences using pure dimensional regularization with d = 4 — 2e and work in the MS 
renormalization scheme. 

The renormalized hybrid-impact-parameter space soft function S{q~ , q~^, b±, fi) is related 
to the bare function Sb{q~ , q~^,bj_) as 

S{q',q^,b_L,fx) = j dui j duJ2 Zg\q' - ui, q^ - uj2, fx) Sb{uJi,uj2,b_L), 

(9) 

where Zs{q~ — Wi, g+ — a;2, A*) is the ultraviolet renormalization constant with an expansion 
around the e — )■ limit given by 

oo ^ 

Zs{uJl,U2,ij) = 5{uJi)S{uJ2) + '^—Zs,kiuJl,UJ2,fJ'), 

k=l ^ 
oo ^ 

Zg\uJi,UJ2, fl) = 5{Ui)S{uj2) + — ^5,fc(^l,^2,/i)- 



k=l 



The Zs k and Zs k can be related to each other from the condition 



(10) 



du[ I duj'2 Zg^{uji — uo[, UO2 — UI2, f^)Zs{<^'i — ^'i-i ^'2 ~ ^21 /^) — ^{'^i ~ uo'l)5{uo2 — Wg)- 

(11) 
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The bare and renormalized strong couplings and as respectively are related by the renor 



malization constant Zg as 



a. 



TT 



2 



(12) 



where = (A^f/12 - llCA/24)/e at NLO. 

The renormalization group evolution of the soft function is determined by the equation 



(13) 



where the anomalous dimension is defined as 



d 



75(a;i,W2,/i) = - I duj[ I du'^Zg {ui - u[,U2 - u}'2, fi) fi—Zs{u}[,U2, fi). 



(14) 



From the finiteness of the anomalous dimension it can be shown that at any order in per- 
turbation theory it is given by 

d - d 
75(^1, W2,/i) = -2as— Z5,i(wi,W2,/i) = 2as— Z5_i(wi,W2,/i)- (15) 

The as expansion of the anomalous dimension is defined as 



75(^1, W2,/i) = ^ 



n=l 



asifi) 



TC 



(n) 



7?^(wi,W2,/i) 



(16) 



Including the as expansion for each pole term in Eq. ([TI1|), we can write an expansion for 

Zs{uJi,U2, fJ,) and Zg^{ui,U2, ^) in ag and e as 



00 00 



Zsiuji,U2,n) = 6{ui)6{u2) + 



1 lasifi) 



k=i j=i 

00 00 



TT 



Z^£^{Ui,U2,fi). 



Z,\uJ^,uJ2,^^) = 5(u;i)5(a;2) + $^5^^[^J'^g(c^i,c^2,/i). 



k=i j=i 



(17) 



The renormalized and bare soft functions have perturbative expansions given by 



Sb{q ,q^,bj 



j=0 

00 J, 

j=0 



Si'\q-,q^,b^). 



(18) 
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The renormalized soft function can be obtained at each order in perturbation theory by 
using Eq. (|I2]) to write the bare couphng in terms of the renormahzed couphng and then 
equating powers of as in Eq. (^. The resulting consistency conditions at LO, NLO, and 
NNLO are given by 



oo -, 
k=l ^ 

S'^'\q-,q+,b^,f,) = 2zg,f,l^Sl'\q-,q+,b^) + f,',^Sl'\q~,q+,b^) 

oo ^ 



+ /io' f dui I du2^\z^gl{q - uji,q+ - uj2, fi) S^'^\ui,U2,b±). 
fc=i ^ 

(19) 

In the MS scheme the renormahzation constants Z^g\, and Z^g \ are determined by requiring 
a cancellation of all pole terms of the RHS above in order to yield a finite result for the 
renormalized soft function. 

A similar analysis can be done for the position-space soft function, and we outline the 
main features below to establish notation. In position space, the renormalized and bare 
soft functions S{b,fi) and Si,{b) respectively are related by a multiplicative renormahzation 
constant so that 

Sib,fi) = z;\b,^)S,{b). (20) 

The anomalous dimension of S{b, /i) is defined as 

d d 
/^^'^(&,/^) = lsib^,b',fi) S{b,fi), Jsib,fi) = -Z5^(6+,6",/i) fi—Zs{b^,b-,fi), 



(21) 



with an expansion in given by 

-a,(;u)i{«) 



oo oo 



75 1 

n=l 

The expansion in as and e of the renormahzation constants are defined as 

1 lasifx)- 

k=i j=i 

1 \asin)- 

k=i j=i 



Zs{b\b-,,) = i + EE^m'^S(^^^-'^)' 
k=i j=i 

z;\b\b-,,) = i+j:j:^[^]'z^^iib\b-,,), 



(23) 



7 



where once again the Zg\{h'^,h~,ij) and Zgl{b'^ ,b~ , fi) coefficients will be related to each 
other by the condition Z~^{b^ , fi)Zs{b'^ , fi) = 1. The perturbative expansion for the 
bare and renormalized position-space soft functions are defined as 



S{b,fj.) 



j=0 



TT 



S,{b) = J2[-^ 

j=0 



(24) 



and the joint expansion in and e of the bare soft function is defined as 



oo oo 



k=l j=l 



1 [a 



(25) 



Following the same procedure as in the case of the hybrid impact-parameter space soft 
function 5(g~, g"*", /i), consistency conditions analogous to Eq. (|T^) can be derived for 
the position-space soft functions S{b,fi) as well. 



B. NLO Soft Function 



We calculate the soft function by inserting a complete set of states using the identity 
l = J2\Xs){Xs\ in Eq. (|6]) to get 



Sib,fi) = ^Y.^r{0\T[SiSn]{b)\X,){X,\ T[45„](0)|0), 



(26) 



and then compute the resulting product of matrix elements in each term above. At LO in 
QCD perturbation theory, only the vacuum state \Xs) = |0) contributes. At NLO, both 
virtual corrections to the vacuum insertion and the single gluon final state \Xs) = 
contribute. In pure dimensional regularization the virtual graphs are scaleless and a non- 
vanishing contribution arises only from the real emission of a gluon into the final state. 

From the results of Ref. |]T9|, the LO and NLO terms of the bare soft function 
<Sb{(l~ , <l~^,bj_) of Eq. dT^ ) are known to be 

Si'\q'-,q-,b^) 



Si'>iq^,q-,bj 



(4vr) 



Ffl 



(g+g ) 



-l-e 



li^l 1 



b]_q~^q~ 



(27) 



We note that g^ are constrained to be positive quantities. Using the NLO consistency 
condition corresponding to the second equation in Eq. ([T9|) and using Eq. (^, the anomalous 
dimension at one loop 7^^^ is given by 



2as^ 


' 1 


' " 


5(g+) + - 


' " 


+ 


TT 


. 




+ 





(28) 
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We note that we have also performed this NLO calculation using an off-shell regulator for 
infrared divergences and have obtained the same results as above. 

In full momentum space, the LO and NLO coefficients of the bare soft function Sb{q,fi) 
are given by 

sr(?) = ^""('!), 

where the + subscript on the delta function denotes that only the positive energy solution 
is taken. It is straightforward to check that the Fourier transform of these coefficients with 
respect to q± give the coefficients of Eq. (pTf ) in hybrid-impact-parameter space. 

Finally, the soft function in full position space is obtained by Fourier transforming the 
result in Eq. (pQ]) to get 

Sl'\b) = 1, 

Sl'\b) = C^^^^e'^^-f^^'^L^ f-e, -e; 1 - e; ^] . (30) 



L = -b+b-fi^e^^'^/A. (31) 



where we have defined 



Upon taking only a time-like component b b^, this is in agreement with Ref. ||5^ (we 
note that /io must be set to unity when making the comparison of our bare results with 
Ref. [|^, as our S^''^ are defined as coefficients of the bare coupling constant rather than 



the renormalized one). The one loop anomalous dimension 7^^'' for the position space soft 
function S{b, jj) is given by 

^f{b,^^) = ^CplnL. (32) 

TT 



III. THE SOFT FUNCTION AT NNLO 



For the calculation of the soft function at NNLO, we again start with Eq. (^) for the 
position-space soft function. Contributions will arise from two-loop virtual corrections to 
the ), one-loop virtual corrections to single-gluon emission (Sj^y), and 

the real emission of two gluons or a quark-anti-quark pair in the final state (5*)^^) so that 
we can write 

S^'\b,fi) = 4'^(6,/i) + 4y(^/^) + sSkb,f^)- (33) 

(2) 

In pure dimensional regulation, the purely virtual contribution Syy is scaleless and vanishes 
and will not be studied in further detail. Representative diagrams for the remaining non- 
vanishing contributions are shown in Fig. |l|, where the dark solid lines represent the soft 
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FIG. 1: Example diagrams contributing the the soft function at next-to-next-to- leading order. 
From the top left diagram and proceeding in a clockwise fashion, the diagrams are respectively 
contributing to the terms S^^^^^,, S^j^CpCA^ ^RV,CfCa ^RkcFN^- 

Wilson lines and the dashed line indicates the insertion of states that puts particles on 
their mass shells. These contributions can be further decomposed according to their color 
structure so that we can write 

S'^^^(6, /i) = Syy{b, /i) + "S*^^^^!^ + S^rr,CfCa + ^RR,CfNp + ^RV,CfCa' (^'^) 

where the subscripts Cp, CfCa, and CpNp denote corresponding the color structures and 
Np is the number of massless fermions. We will consider each of these contributions sepa- 
rately, presenting both the final result and the relevant technical details. We will first present 
results for the bare soft function and then discuss renormalization and the the extraction 
of the anomalous dimension. As checks of our results we will use both the comparison with 
the 6 — )• 6" limit of Ref. and the constraint of non-abelian exponentiation [^, We 



will present results in both position-space and in the hybrid-impact-parameter space. 
A. The bare soft function 

We begin by first presenting the results for the various contributions to the bare NNLO 
soft function in position space. The NNLO calculations are first performed in momentum 
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space and then converted to position space by inverting Eq. (^. The most difficult integrals 
occur in the double real-radiation graphs; the necessary integrals are presented in Appendix 
^ The results for these integrals, which typically contain hypergeometric functions, must be 
Fourier transformed to position space. To do so, we ffist series expand the hypergeometric 
functions in their argument, perform the integrals for each term of the series, and then 
identify the resulting sum. The relevant integrals for this procedure are presented in the 
Appendix. 

The results for the various non- vanishing contributions to the bare soft function Sb{b,fi), 
in pure dimensional regularization are 



(2) 

,RR,CpNp 



^(2) 
bfRFt^C p 



(6) = -C^n /'^^.^ ^—^ ^e-2^^>o''^'Si^i(-2e,-2e;l-2e;-^ 

^ ^ gS 8(1 -2e)(3-2e) in,, , 



4e4 



2e7i5 ,,-4£ r2e 



3F2(-2e, -2e, -2e; 1 - e, 1 - 3e; ^ 



+ 2i^i(-2e,-2e;l -3e 

^2 



■'b,RR,CFCA 



+ 2 



16e4 
r2(l-2e) 



6+6- 



+ 0(6) 



^PL.„..(6) = C^C^^4^e"^-^>n ^-L^- | 7?!! ~ . 2^1 (-26, -2e; 1 - 2e; 



;i- 26)(3 - 26) 



6+6- 



r(i -36)r(i 



3F2(-6, -6, -6; 1 - 6, 1 - 36; 1) 



X 3F2(-26,-26;l-26;l-6,l-36;^)- 2 2Fi(-26,-26;l-36 ^ 



6+6- 



' 6+6- 



+ 0(6), 

= —CaCf- 



r'(i-^)„-2, 



^_2e7s^-4e^2e ^p^^ _ 2e)r2(l + e) C0s(7r6)] 



X 2i^i (-26,-26;! 



6+6- 



(35) 



For simplicity of presentation, we have expanded some of the results above only up to 
the needed order in 6. It is straightforward to conffim that in the limit 6_l — )■ 0, these 
expressions are identical to those obtained in Ref. |]5^. In addition, the full position-space 
result is constrained by non-abelian exponentiation |5^, which requires the terms at 
NLO and NNLO to obey the relation Sf^^^ (6) = [^^^^6)] ^ /2. This relation can be checked 



B3, |BTO]). These two 



through the finite order in e using the expansions given in Eqs. 
checks are strong indications of the correctness of our results. 

We now present the results for the soft function in the hybrid impact-parameter space, 
that appears directly in the formalism of Refs. 



These results can be obtained 



through the Fourier transform of the position-space results. As in the case of the position- 
space soft function, we separate the results into several components dictated by the color 
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and cut structure, and obtain the following expressions: 



X oi^i(l-2e 



er2(l-2e)2(l-2e)(3-2e) 



o(2) 

'^b,RR,CFCA 



i2e 





4 


e2r2(l 


-2e) 


-e,l- 


-3e;- 







), 

b±'^q~^q 



oi^i(l-3e; 



6_L g+g" 



+ 0(6), 



X oi^i(l-2e 



e2r2(l-2e) 



{q-^q- 



+ 



/ (2-e)(3-e) 
4(1 -2e)(3-2e) 
r2(l -2e) 



2r 1 - 3e r 1 - e 



^ 2 + - 

X 3F2(-e, -e, -e; 1 - e, 1 - 3e; 1) 1^2(1 - 2e; 1 - e, 1 - 3e; - ^ ) 



4 



-- oFi(l -3e; ^ 



0{e) 



0(2) 



X oFi(l-e;- 



r2(i_e) r(l-2e)r2(l + e)cos(7re), + 



r2(l-2e) 

6_L q^q . 



2e2 



-(g+g" 



4 



(36) 



The C*!^ piece again satisfies non-abelian exponentiation, but with a convolution in momenta 
rather than a simple product as explained in the next section. 

B. Renormalization, exponentiation, and the finite soft function 



The soft function is constrained by non-abelian exponentiation Q so that one can 
write the position-space soft function as 



S{h,^^) = exp{.(6,^)} = exp | ^ (^)%(")(6, /i) 



(37) 



with appropriately defined s{b, fx) and s*^"^(6, /i). The bare and renormalized position-space 
soft functions are related as in Eq. (^0]). In exponentiated form, this relationship can be 
re-expressed as 



S{b,iJ,) = exp{z5(6+,6 , fi) + Sb{b)} 



(38) 
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where 



oo 2n 



n=l m=l 



7r 



-("•) 



oo 2n 



Sb 



n=l m=— oo 



vr 



(39) 



s['^(6). We have introduced a similar exponentiated 



with = ^ , and 

form for the renormahzation constant, 

= exp {^5(6+, 6", /i)} . 



(40) 



The bare and renormahzed strong couphngs are related by Eq. i\i2\) . The non-abelian ex 
ponentiation theorem implies consistency relationships between the coefficients in Eqs. (p3l 



and ( |25| ) and the coefficients and Zg^. 

The renormalization-scale independence of the bare soft function determines the anoma- 
lous dimension in terms of zs as 



An) 



fi-^S{b,fi)='ys{b'^,b ,fx)S{b,n), 
an 



75(6+, 6 = n^zs{b'^,b ,/i). 
an 



(41) 



The finiteness of the renormahzed soft function S{b, fi) determines the renormahzation con- 
stants in Eqs. (|38[) , (0). The LO and NLO anomalous dimension contributions are given 
by 



7s 

(2) 

7s 



CfNf 



5, ^ 14 

— InL \ 

9 27 36 



+ CfCa 



\18 6 



101 llTT 



27 



72 



(42) 



By comparing the coefficients of the expansion of the RHS of Eqs. ( pT)) and (^) after 
renormahzation, the NLO and NNLO terms in the exponent of the renormahzed soft function 
S{b, n) in Eq. (^) are determined to be 



.«(6,/x) 
.(2)(6,/i) 



- n^LH hLi2(— 4- 

2 12 ^^6+6- 



CpATp <^ In^ L In^ L - — + -Li2( 



Lio 

18 



36 
b^ 



36 



7 L 



27 6 



'b+b- 
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CfCa <^ — In^ L + In" L + 



72 



67 TT^ 



72 24 

4 



101 7 



InL 



41 



162 432 



Svr^ ^ C(3) 



36 



54 



60 7 677r2 11 vr 
^324 ^ "86? ~ 72^*^ '* ~ 48 



67 7r2\ . ^ bl ^ 11 . ^ 6i 
36 12/ ^^6+6-^ 12 ^^b+b' 



InL 

1. 



:Li4( 
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-lis,.(A,)_s„(A,)-lLi5,i,)}. ,43, 

Here, Li„(x) is the usual polylogarithmic function, S„,p(x) is Nielsen's generalized polyloga- 
ritlim, and we remind the reader that L = — 6+6^yU^e^'^^/4. 

Non-abelian exponentiation also holds in the hybrid-impact parameter space but with the 
multiplication operation replaced by a convolution in the light-cone momentum coordinates 
so that 



S{q ,q^,b^,ij,) = S^°\q , g+, (g) exp ,g+,6_L,/i)} 
where the perturbative expansion of the exponent is defined as 



(44) 



n=l 



(45) 



It should be understood that the first term in the expansion of the exponential is 6{q )6{q^). 
The first few terms obtained after expanding the exponential in Eq. (|4^) take the form 



5(g-,g+,6x,/x) = 5(°)+5(°)® [s + ^ s ® s + ■ ■ ■ ] , (46) 

where we have suppressed the arguments of the functions on the RHS. The convolution 
product (X) in the light-cone momentum coordinates is defined as 

fiq~,q^)®9iq~,q^) = j dui j duj2 f{q' - uji,q^ - UJ2) g{uJi,uj2)- (47) 

Similarly, the relationship between the renormalized and bare soft function in exponentiated 
form also involves convolution products and is given by 



S{q ,g+,6_L,/i) = Sy{q ,q^ ,b±) exp {zs{q ,q^,l2) + Sb{q ,q^,b±)}, 
where the expansion in as for the terms in the exponent are defined as 

00 2n ^ ^{n) 00 2n _M 

^5(g~,g+,/^) = $^$^(^) s,{q-,q+,b^) = Y^ Yl ( 



(48) 



n=l m=l n=l m=—oo 



TC / e' 



(49) 



with z^gl^ = z^gl^{q ,g"^,/i) and s^^^ = s^f^l^{q ,q~^,b±). The anomalous dimension and the 
renormalization constants again satisfy 

lji-^S{q~,q^,b±, n) = J dui J duj2^s{.q~ - ^i,q'^ - ^2, jj) S{uJi,uj2,b^, jj), 

14 



(50) 



For simplicity of notation, in the rest of this section we use the following definitions 
T{z) = Jo(.)(ln(|)+7i.)-flo(.), 



-^0,0 



6{q-)6{q' 



1 



+ 



1 



2 
/i 



1 


"ln(g+//i)" 




fi 


_ q+ 


/fi 




+ 


1 






fi 


. q' 


-/fi 




1 


' ' 








[q^l 


+ 


.q 





siq-) + 



fx 



6{q-) + 



fi 



Kq^). 

" ln(g-//i) 

. q~/^^ 

^InVM 

q-/fi 



ln(g+//i) 
q+/fi 



q~ 



/i2 



Mq //^) 
q~/^^ 



ii 



(51) 



We note that T(0) = 0. The LO and NLO anomalous dimension contributions in the 
hybrid-impact-parameter space are given by 



is\q ,q^,f^) 
is\q~^q^^f^) 
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<^ C\N^ 


[-( 












27 


36 


CfCa 








llTT^ 


( 






72 






J7 





1 -^0,0 + 



+ 



/67_7r2 
Vl8 6 



(52) 



The same results for the anomalous dimension can be derived using Eqs. (JT^) and (|T9|) 
without using the exponentiated form of the soft function, which provides a consistency 
check on our calculation. The results for the NLO and NNLO renormalized coefficients of 
Eqs. (H3) and (E3p are 



s(^)(g+,g ,6±,/i) 
s^^\q^,q~,b±,fi) 



CfNf 
5 



7T 

12 



^0,0 + ^0,2 + Mb±Vq^q')Li,i^ 



7 



TT 



Lo,2 + 



" ~ i^^^O ^ I 27 " 36 ' 
\T{h^^Jq+q- 



0,3 



6 



18 



-''o(^±\/ 



^1,1 



+^Jo(&±Vg+g-)^o,2| + 



15 



(53) 



where Jn(x) is the standard Bessel function. 



IV. CONCLUSIONS 



We have described a computation of the exclusive soft function for Dreli-Yan production 
of electroweak gauge bosons through next-to-next-to-leading order in perturbation theory. 
This object is required for the resummation of low pt logarithms through next-to-next- 
to-leading logarithmic accuracy in the SCET approach of Refs. [|I^, |TP|. Results for both 
the anomalous dimension and the finite soft function have been presented, and all relevant 
technical details have been explained. Adapting these techniques to the computation of 
the soft function that appears for gluon-initiated production of a Higgs boson should be 
straightforward. We expect that the exclusive soft function will have further applications 
in precision studies of differential distributions within SCET, and that our result will be an 
important step toward enabling these future studies. 

We conclude with a few comments on the SCET approach with unintegrated distribution 
functions to the low pt distribution, in which the need for this soft function first arose. 
In the standard TMDPDF approach [p5|-|29|, |3T1^5[, rapidity divergences arise in pertur- 



bative computations that require additional regulators beyond the standard dimensional 
regularization. The need for similar regulators also arises in SCET approaches to related 



observables [0. In our approach, the Impact-parameter Beam Functions (iBFs) and the 
Inverse Soft function (iSF) are more differential in momentum coordinates than the corre- 
sponding objects in the TMDPDF formalism. As a result, in perturbative computations of 
the iBFs and the iSF the rapidity divergences are regulated by the physical kinematics of the 
process. An investigation of the relationship between the our approach with iBFs and the 
iSF and the TMDPDF formalism is worth pursuing in future studies. However, some recent 



results [£8| based on the TMDPDF formalism do not include the soft function needed for 
a proper treatment of soft radiation and lack operator definitions, preventing any rigorous 
field-theoretic interpretation and rendering comparison with our results difficult. 

We have performed the first step needed for NNLO studies of low pj^-transverse momen- 
tum distributions, using the SCET approach with unintegrated distributions functions, by 
computing the relevant exclusive soft function at this order. These results are also the first 
step towards achieving low-pr resummation at NNLL accuracy. We look forward to the 
further development of our approach to low-px distributions. 
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Appendix A: Auxiliary integrals 

We compile here several integrals that appear during the course of our calculation. We 
define 

[dk] = d% d% S+{kl) S+{kl) S'^iq - fci - k2). (Al) 
It is straightforward to derive the following by direct integration: 

The following integrals are useful in performing Fourier transformations: 



I 



I 



I 
I 
I 



^^i-e ^,+,-u-.+(m+n)/2 r(l + m/2)r(l - 6 + n/2) 
^ ' r(l-e)r(2-e + (m + n)/2) 

x.F.(l-e+^;l-e,2-e+^;-^), (A6) 

^d-2q^ ^-iq^-b^ q-2-4e ^^^^2^ 2Fi(-e, 1 - e; ^ 



= TT-^ oFi(l - 36; -^^^) + O(e^), (A7) 

e-m-6x q-2e ^-2-2. ^^^2) ^ _ ^. ^) 

= TT^- (^^^")"'^Z^T^rT^ 1 - 6, 1 - 36; -^^^) + 0(6^), (A8) 

d'^-V e-'^-'^-g-2-2^gl2^M?')2Fi(-e,-e;l-6; ^ 



(^+^-)-3e^il_Jl_ ^p^^_^^ 1 _ e, 1 - 3e; 1) 
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X 1^2(1 - 2e; 1 - e, 1 - 3e; 



4 



0(. 



(A9) 



+ — ?,2 



, ai, . . . , a^; bi, 



b+b- 



(AlO) 



The + subscript on the delta and step function denotes that only the positive energy solution 
is taken. Although some integrals are only correct up to O(e^), they become exact in the 
limit of 6x — )■ or in the presence of S{q^). 

Appendix B: Hypergeometric expansions 



For completeness, we present here several expansions of hypergeometric functions that 
we found useful in our analysis. These can be simply obtained using the series expansion 
of the hypergeometric function, expanding the resulting Gamma functions in e, and using 



known techniques for summing the resulting series in terms of known functions |59 



pFq{ai, a2, . . . , ttp] bi, 62, • • • , bq] 0) 
2Fi(a,6; c; 1) 



oFi(l - ae; -—] 



iF2(l-2e;l-e,l-3e;--) 



r(c)r(c-a-6) 
T{c-a)T{c-b) 
17 



,F2i-e, -e, -e; 1 - e, 1 - 3e; 1) = 1 - C(3)e^ 



360 



r(l - ae) { Mz) + ae Uz) ln(|) - ^Yo{z) 



+0(6^)} 

r(l-e)r(l-36) 

r(l-2e) 
+0(6^)} 



(Bl) 
(B2) 

(B3) 
(B4) 



Uz) + 2e 



Jo(^)ln(-)--ro(^) 



(B5) 



iFA—ae, —ae] 1 — ae; z) 



iFi(-2e, -2e;l-3e;z) 



1 + ah^hhiz) - ah^ [81,2(2) - Lh{z)] + 
aV [Li4(z) + Si,3(^) - S2,2(^)] + 0(e5) (B6) 
1 + 4e2Li2(z) - 4e3 [81,2(2) - 3Li3(z)] + 
2e^ [Liliz) + 18Li4(z) + 281,3(2) - 1082,2(2)] + 
O(e^) (B7) 
,Fi(-2e, -2e; 1 - e; z) = 1 + 4e2Li2(2) - Ae^ [381,2(2) - 113(2)] + 

2e^ [Liliz) + 2114(2) + 1881,3(2) - 1082,2(2)] + 
0{e') (B8) 
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2Fi(-2e, -2e, 1 - 2e; 1 - e, 1 - 3e; z) = 1 + Ae^U2{z) - 8e^ [Si,2(^) - Us{z)] + 

[Liliz) + 10LU{z) + 8Si,3(-2) - 12S2,2(^)] + 
0(6^) (B9) 
3F2(-2e, -2e, -2e; 1 - e, 1 - 3e; ^) = 1 - 8e=^Li3(z) - 16e^ [2Li4(^) - S2,2(;2)] + O(e^) 

(BIO) 

Jo{z) and Yo{z) are the standard Bessel functions, and Li„ denotes the standard polylog- 
arithmic functions. In addition, Nielsen's generaUzed polylogarithms, denoted by Sn,p{z), 
appear. In the final finite results for the soft function, only two of these functions appear. 
They can be exchanged for the standard polylogarithms using the following identities: 

S,,2(2) = + ln(l - z)U2(l - 2) - Li3(l - 2) + C(3), (Bll) 

4- (^'2' 
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